MAN/HEP/2006/29 

Hybrid Inflation Exit through Tunneling 



o 



Bjorn Garbrecht 
School of Physics & Astronomy, The University of Manchester, 
. Oxford Road, Manchester Ml 3 9PL, United Kingdom 

O 

Thomas Konstandin 
Department of Theoretical Physics, Royal Institute of Technology (KTH), 
AlbaNova University Center, Roslagstullsbacken 21, 
106 91 Stockholm, Sweden 

CM 

Abstract 

For hybrid inflationary potentials, we derive the tunneling rate from field config- 
CM ■ urations along the flat direction towards the waterfall regime. This process competes 

with the classically rolling evolution of the scalar fields and needs to be strongly sub- 
dominant for phenomenologically viable models. Tunneling may exclude models with 
a mass scale below 10 12 GeV, but can be suppressed by small values of the coupling 
constants. We find that tunneling is negligible for those models, which do not re- 
'"q! quire fine tuning in order to cancel radiative corrections, in particular for GUT-scale 

i 1 SUSY inflation. In contrast, electroweak scale hybrid inflation is not viable, unless 

the inflaton-waterfall field coupling is smaller than approximately 10 -11 . 

> 
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The slow roll paradigm of inflation [1] requires the scalar potential to be flat to such an 
extent, that the Hubble expansion causes an overdamping of the evolution of the inflaton 
field. This has the consequence, that the kinetic energy of the inflaton is negligible, and the 
equation of state of the dominant component of the Universe is approximately the same 
as for vacuum energy. When realizing slow roll inflation within single field models, one 
encounters the problem of reconciling the flatness of the potential, its comparably large 
magnitude and the wish to keep the vacuum expectation value (VEV) of the inflaton below 
the Planck scale. 

A possibility to address this problem is the hybrid inflation |2J mechanism, where the 
slowly rolling inflaton triggers a waterfall field to rapidly roll down the potential and to 
terminate inflation at some critical point. The direction along which the potential in- 
creases towards large values driving inflation and the direction of the slow-roll are therefore 
separated. 

When comparing different models of hybrid inflation at the same scale, that is with 
the same value of the potential, it is clear that in a model which has a flatter direction 
for the inflaton, a certain comoving scale leaves the horizon when the inflaton is closer to 
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the critical point than in a model with a steeper direction. Imagining the limiting case 
of a completely flat direction, the classical field dynamics suggest that inflation may last 
infinitely long with the inflaton being arbitrarily close to the critical point. However, within 
quantum theory, metastable configurations eventually always decay to the one of lowest 
energy. We therefore expect that in hybrid inflation, a field configuration along the flat 
direction may tunnel to form a bubble containing a field configuration in which inflation 
ends and the scalar fields rapidly assume the true vacuum state. It is the purpose of this 
study, to estimate this decay rate, compare it to the classical field evolution and to specify 
for which model parameters tunneling is a non-negligible effect. 



2 Tunneling during Inflation 

2.1 Tunneling without Barriers 

The semiclassical theory of tunneling for scalar field theory is developed by Coleman and 
Callan in [3J 0] . We consider the Lagrangian 

c = \{d&W<p)-v{<p), (i) 

where the field <p is initially located everywhere in space at the point ip + , which corre- 
sponds to a false vacuum or a classically metastable configuration, and where we normalize 
V(ip + ) =0. In order to calculate the decay rate, one proceeds by solving the classical 
Euclidean equation of motion 

where the prime denotes a derivative with respect to if. We assume that the solution takes a 
spherical symmetric form in Euclidean space and write g = \x\. In order to understand the 
properties of the solutions to this equation, it is most useful to recall that it corresponds to 
the equation of motion for a one-dimensional particle moving in the potential V(<j>) turned 
upside down and with a friction term (3/ g)(dip/dg), which implies infinite damping at 
g — and vanishing damping when g — > oo. 

The instanton solution, which obeys the boundary condition ip(oo) = (p+ is called the 
bounce, and we denote it by Tp(g). It uniquely determines a release point ip T , at which 
V(ipi) < V((p+) and which satisfies dlp/dg — at Tp = <f r and g = 0. Physically, y? r is the 
initial value of the scalar field inside a nucleating bubble, from which it starts to evolve 
classically. 

Having found the bounce solution, we can compute its Euclidean action 



S E = 2vr 2 / g 3 dg 



1 fdlp x " 



2 V dg 



(3) 



which is used to obtain the tunneling rate T per volume V as 

r = si ( det'[-d" + v"m} y 1/2 c _ SF 



V 4tt 2 Vdet[-<9 2 + V"(<p+)]J 



e~ bE , (4) 
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where the prime at the determinant indicates the omission of the zero eigenvalues. The 
evaluation of the determinants is a quite costly task, and we follow the common prac- 
tice jSl EH Hj to estimate their values from the parameters of the particular theory under 
consideration. Indeed, the results we present justify this procedure a posteriori. 

We intend to apply this theory of tunneling to hybrid inflation, which is implemented 
by the generic potential j2] 

V(a, 0) = V (a, 0) + V L (a) = - — 2 + — + -g 2 0V + V L (a) . (5) 

This potential is almost flat with respect to the inflaton a along the direction where = 0. 
The flat direction is lifted by the contribution Vl(ct), where we normalize Vl(0) = 0, which 
causes a to classically roll down the potential from larger to smaller values. Inflation ends 
shortly after a reaches the critical value 

a c = —. (6) 
9 

At this point, the mass square for the field changes its sign from positive to negative and 
the inflationary valley turns into a ridge. The field then quickly evolves away from zero 
and the fields eventually assume the values 

Tfl 

°o = , 00 = -j= , (7) 

where V(ao, 0o) = and inflation is terminated. Due to the transition from valley to ridge, 
from which the fields fall, this is called the waterfall mechanism, and we denote the area 
where a < a c as the waterfall region. 

Returning to the question of tunneling, we note that the hybrid potential © does not 
have any local minima but the global one (|7jl. Therefore, there are no false vacuum configu- 
rations possible, and it may appear that the theory of tunneling and bubble nucleation does 
not play any role for hybrid scenarios. However, as already mentioned in the introduction, 
we can imagine the case Vl(ct) = and wonder whether a configuration with a > a c is 
stable. 

Quite similar situations are discussed by Weinberg and Lee j^j, and they point out that 
a bounce solution can exist in some cases without a potential barrier between the initial 
point and the global minimum of the potential 1 . The necessary condition for the existence 
of a bounce is not the presence of a potential barrier, but of an energy barrier, constituted 
by the potential and a contribution from the gradient terms of the bubble wall. Therefore, 
a false vacuum is not required to exist for tunneling to be a relevant process. A very 
instructive example is given by the potential 



V{<j>) 



0, < 
-k<f> > 



1 Linde has given an earlier example of upside-down <j> -theory, where tunneling can occur [H]. See also [7] 
for a more recent related discussion. 
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Classically, if the field is positioned on the plateau at the position = — A0, the system 
would be stable, while quantum-mechanically, it turns out to be unstable due to tunneling. 
The existence of a corresponding bounce solution can be understood from the Euclidean 
equation of motion If the field is released at rest when g = and = r > 0, it will 
accelerate in the upside-down potential until = and then asymptotically come to rest 
again at = —A(f) due to the damping term. This bounce solution therefore describes 
tunneling from the metastable position = — A0 on the plateau to nucleate a bubble with 
the vacuum expectation value = r inside. The rate for this to happen is calculated to 
be [SI 

r 4 2A . / 32vr 2 A0 3 



v V^ ex H-^^J' (9) 

where C is a constant of order one, which can in principle be determined by evaluating the 
determinants in Eq. (0J). This result is apparently already very useful in order to estimate 
whether for a given inflationary model, it is in order to worry about tunneling. If the cube 
of the distance from the region where inflation takes place to some other point of lower 
potential is of the same order or smaller than the derivative of the potential at that point, 
the bounce action can be of order one and tunneling sizeable. Similar to this example, for 
the hybrid potential Vo, Eq. (j5J), bounce solutions exist that start in the waterfall region 
and come at rest on the flat direction where = and a > o c . 

One may argue that the potential during inflation is not exactly flat and that therefore 
the formula Q for the tunneling rate does not apply. We follow however the argument 
of Weinberg and Lee, that taking the motion of the inflaton field or the lifting of the flat 
direction into account will only reduce the action of the tunneling process. For calculating 
the bubble nucleation rate in the hybrid model, we therefore determine the bounce solution 
for the potential Vq and neglect the effect of Vl. This way, we obtain a lower bound for 
the tunneling probability, which still allows to derive constraints on the parameter space 
for hybrid inflation. 

2.2 Numerical Results 

We now determine the bounce action for the hybrid potential (0) as a function of the 
distance of the inflaton from the critical point, 

Act = a - a c . (10) 

While it is not possible to find analytic bounce solutions, one can reduce the problem con- 
siderably by making use of the scaling properties of the potential. Inspecting the Euclidean 
equations of motion (J2J) for the hybrid case, 

d 2 a 3 da 



dg 2 g dg 
d 2 6 3, 



2 j.2 

9 o 



()i>2 + = ~ m( t> + r^<i> + 4A0 d , (n) 

we see that they are left invariant under the following rescaling: 

A -> Xk, pn~ 1/2 , m m K 1/2 , g -> gn 1/2 . (12) 
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The bounce action (J3J) then transforms as 

S E (X,m,g,Aa/a c ) = A^S^l, m/VX, g/VX, Act/ct c ) . (13) 
Another rescaling leaving the equations of motion (|TT| invariant is 

m — > urn, a — > k<t, — > K(f), p — *> kT 1 p . (14) 
This reveals that 5*£;(A,m, g, Aa/a c ) does not depend on m, 

S E (X,m,g,Aa/a c ) = A^S^l, m , g/VX, Act/ct c ) =: \~^x{g/V\~, Act/ct c ) , (15) 

where m is arbitrary. 

We now determine the function \ numerically. In general, finding bounce solutions 
can be very complicated for multi-dimensional problems, or at least time consuming. Two 
algorithms, that can be applied to a wide range of problems, have been presented, e.g. 
in Refs. [HJ E]- These algorithms are not immediately applicable to our problem, since 
they have been designed for the case of tunneling with potential barriers. Fortunately, for 
two-dimensional problems, one can resort to scan procedures, which we apply here. First, 
we fix the starting point of the configuration (cto, 0o) and solve the equations of motion 
by integration. For late times, the solution can behave in two qualitatively different ways. 
The first possibility is that a always stays smaller than ct c , and oscillates around zero. 
In this case Cto was chosen too small. In the second case, a is finally larger than ct c and 
the upside-down potential is hence unstable in the 0-direction. Depending on the initial 
point, the configuration then behaves usually as <fi — > ±oo, when p — > oo. These two cases 
correspond to the 'over-/undershooting' of the one-dimensional problem. Keeping 0o fixed, 
while varying o"o using the 'over-/undershooting' method, leads thus to a bounce solution. 

In Fig. [H we plot the function ^(g/v'A, Aa/a c ), obtained by the above procedure, for 
different values of g/VX. The results show for small Act -C mg/X a scaling according to 
g~ 2 Aa, which we explain below. The numerical coefficient turns out to be 

S E w 158 x —\. (16) 

To ensure this scaling behaviour even for very small values for the coupling constants, we 
will give an analytical upper bound for the Euclidean action in the following. With above 
insights, we can proceed with further simplifications of the problem. Sizeable tunneling 
may only occur when the inflaton a is close to its critical value, cf. Fig. ^ Therefore, we 
assume 

Act < ct c . (17) 

In order to obtain a lower bound on the tunneling rate, we impose the instanton to follow 
a straight trajectory in (ct, <ft) space. The exact solution along a curved trajectory has a 
lower Euclidean action and therefore corresponds to a larger tunneling rate. The trajectory 
is parameterized by 

(f> = aw , 

ct = Act + — - Vl -a 2 w, (18) 
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Figure 1: The two-dimensional numerical result for the function xid/v^i Aa/a c ) for the values 
g/VJ = V2 (dashed), g/V\ = 1 (dot-dashed) and gj\f\ = l/y/2 (dotted). 



where a G [0; 1] is a free parameter that will be determined by minimizing the action. 
Along this trajectory, the potential (jSJ) close to the critical point takes the form 



1 m 



V = -— - aVl - a 2 gmw 3 + a 2 gmAaw 2 + O (w 2 Aa 2 , w 4 ) . (19) 
4 A 

We now determine the value of the parameter a, for which the Euclidean bounce action is 
minimal. For that purpose, we consider the potential 

V = -aw 3 + (3w 2 . (20) 

By rescaling arguments, one obtains that the corresponding action has to scale as 

(21) 



(3 Aa 1 
a 2 mg a 2 (l — a 2 ) 



and is minimized for a = 1/ y/2. This explains the scaling behaviour for small Aa observed 
in ((TBI) . The comparison with Eq. (fTH|) yields for the linearized case 



Se ~ 

x(g/VX,Aa/a c ) ~ 



Act 1 



(22) 
(23) 



For the choice a = l/y/2, neglecting the w 4 terms in the approximated linearized poten- 
tial (fTUj) is justified when 



c A ' 



(24) 
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and the w 2 Aa 2 terms are subdominant if 



Act < 2cr c . (25) 
Numerically, we find for the constant of proportionality 

S E = 182 x — 1 , (26) 

c c g 2 

where the larger factor of proportionality when compared with (|T6*j) is due to the fact that 
we are restricted to the linear path and therefore miss the minimum of the Euclidean action 
in the two-dimensional field space. We also note that the point w r , from which the field w 
in the bounce solution is released, scales according to 

Act m , 

w r = 8.2x = 8.2 x Act. (27) 

o" c 9 

Notice that a small Euclidean action, Se < 1, automatically ensures the requirements in 
Eqs. ()24)1 and ()25)1 and hence the validity of the approximation in Eq. (|T9|) . if g, A < 1. 

Finally, when assuming m to be of order of the Grand Unified Scale 10 16 GeV or less, 
all scales in the problem are larger than the Hubble rate 2 



8ttV 



3m P1 



where mpi = 1.22 x 10 19 GeV denotes the Planck mass, such that gravitational effects can 
be neglected jTU] . 



3 Bounds on Specific Models 



We estimate the relevant values for Act using the standard slow-roll dynamics of the inflaton. 
When the expectation value of the inflaton, at a certain instant during inflation, takes the 
value ct = cr e , the number of e-foldings N e that will elapse until inflation ends is calculated 

as 



N P 



Hdt 



8tt 



V 



m 



pi 



da 



dV/da ' 



(29) 



where we have used the slow-roll approximation 3Hda/dt = —dV/da. One important ob- 
servational constraint is the amplitude \fP-ji of the power spectrum of scalar perturbations 
for the scale fc, that exits the horizon when ct = ct p , 



k 16 V 3 / 2 
6 m P1 dVj da 



(30) 



2 The displacement Act exceeds the Hubble rate as a consequence of imposing the small observed 
value (|3T)> on the the amplitude of the scalar perturbations (|3U)l . 
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Here, we impose the normalization [TT] 



= 4.5 x 1(T 5 (31) 
at k — 0.05 Mpc" 1 . This scale exits the horizon at 

* = 50 4 log '«io4v + 5 lo& »io^v' (32) 

Since k = 0.05 Mpc -1 corresponds to multipole moments around I = 700, the largest 
angular observable scales have exited the horizon about six to seven e-folds earlier. 

A very conservative estimate for Act and therefore the tunneling rate is therefore ob- 
tained by setting N e = 60 and 

Act = ct c — ct c . (33) 

We use this value to compute the Euclidean action (|2l)j) and to estimate the tunneling 
rate (@J). The latter is to be compared with the expansion rate during inflation H, e.g. the 
number of non- inflationary bubbles nucleated per expansion time in one horizon is given by 
T/(VH A ) and should be much less than one. An interesting, but difficult question would be 
to quantify how much less. Due to the exponentially strong dependence of the tunneling 
rate on the model parameters, we omit a discussion of this question by the same token on 
which we do not evaluate the determinants in Eq. (@J). 

We furthermore remark that it appears very likely that for viable inflationary models, 
one has to impose that tunneling also does not occur at much lower values of N e than 60. 
The nucleation of non-inflationary bubbles would lead to very large density perturbations 
on small scales, which induce the production of primordial black holes , which is strongly 
constrained observationally ^3]. We do not discuss this possibility here any further and 
just explore the conservative bound. 

3.1 Blue Model — Quadratically Lifted Flat Direction 

In the seminal work j2] , hybrid inflation is implemented by a quadratically lifted flat direc- 
tion, through the effective potential 

V h {a) = \mla\ (34) 

Due to the positive curvature of the potential along the flat direction, the scalar pertur- 
bations are predicted to be blue tilted, which is characterized by a scalar spectral index 
n s > 1. Using (J2T?|) and the basic potential (jSJ), we can solve for 

CT e = -exp<^- , (35) 



g \2ir m 4 
while the amplitude of the power spectrum (j5U|) is given by 

'2 gm 5 



3 71 A 3 / 2 TOp 1 m^ ' 
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The latter two equations can be solved for m a and a e by assuming that the exponent 
in (J35J) is small, approximating a e ~ cr c , and justifying this a posteriori. We find 



2 9 I 2vr m 5 



and 

fi- ^ { y/QirXPn m Pl 

such that 

1 m 2 



m q m . 
a c = - exp <^ 7== A^ e ^ , (38) 



Aa « / _ — N e (39) 

V6vrAP-^ Topi 

Inserting these into (J2*H|) and using (pTT]) yields 

42 m Ar 9.3 x 10 5 m Ar 

S E = — ; N e = -= N e . (40) 

gVXPn mpi g^/\ m P i 

We now discuss the self-consistency of the above results. For the approximation of the 
potential V by expression ()19)) to be valid for the bounce solution, we have to fulfill the 
relation (|2*4*j) with w = w T . Using (|2*7j) and (|39*j) with N e — 60, we find the bound 

m«2.2xl0^m P i. (41) 
v A 

This condition also ensures the validity of the assumption Act cr c , in particular that the 
exponent in (|H8jl is much smaller than one. 

In order to summarize these results, we present Fig. 121 A reasonable estimate of the 
tunneling rate is given by 

y = iJ» 4mVSE ' (42 > 

since gm is the smallest dimensionful scale occurring in the approximate potential ()19j) . 
We compare the decay rate with the Hubble rate (|28jl . since (r/V)/i/ 4 is the number of 
bubbles nucleating in one Hubble time within a Hubble volume. Note that for the range 
of m for the individual graphs of Y/V in Fig. the consistency condition (|41j) is met. The 
wide range of orders of magnitude covered relativizes the importance of the prefactor of 
the exponential in the expression for T/V, in particular the determinants in Eq. (|4"|). Also 
the precise bound on the tunneling rate loses importance due to its strong dependence 
on m after it has reached is maximum. As a conservative requirement, we may impose 
{T/V) < H A . A bound which is stronger by a few orders of magnitude might be in order 
to accord with observation, but has no significant impact on the tunneling bound on m. 

An important implication to be read off from Fig. (J2J is that for the TeV-scale as a 
special scale of interest, g has to be smaller than at most 10 -11 in order to avoid a fast end 
of inflation through tunneling, provided A is of order one. Besides by suppressing g, we see 
from Eq. (|40|). that also small values of A serve to suppress the tunneling rate. Choosing 
this option however leads to expectation values for after inflation and a during inflation, 
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Figure 2: Hubble expansion vs. tunneling decay in the blue model. The plot shows log 10 R f° r 
R = H 4 (solid) and R = T/V for g = 0.5 (dashed), g = 5 x 1(T 6 (dot-dashed) and g = 5 x 1(T 12 
(dotted). We have chosen A = 0.5. 



which are much larger than m. If one considers m as a cutoff scale or to be closely related 
to a cutoff scale of an effective theory, this is undesirable. 

As a curiosity, we note that we rule out a particular choice of parameters used as an 
example in the original work on hybrid inflation [2], m = 1.3 x 10 11 GeV, g 2 = A = 0.1. 
In this case, T/V = 1.3 x 10 33 GeV 4 , whereas H A = 1.7 x 10 14 GeV 4 , indicating that 
T/(Vif 4 ) = 7.6 x 10 18 non- inflationary bubbles are nucleated during one expansion time 
within a horizon. 



3.2 Red Model 

Since the WMAP3 data strongly prefers a red-tilted scalar spectral index n s , with the best- 
fit value given by n s ps 0.95 |TT], we also study models with a negative curvature along the 
flat direction. A simple possible realization of these is given by 

VL = A 3 a - ^m 2 a a 2 . (43) 

During inflation, the inflaton takes values in between a c and the maximum of VL, which is 
located at a = A^/m 2 . This translates into the requirement 

771 A 3 /AA\ 

-<—■ 44 

g m i 

Note that this model has an additional parameter when compared with the quadrati- 
cally lifted model, which is fixed by imposing the value of the spectral index of the scalar 
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perturbations n s = 0.95 as an additional constraint. It is calculated through the slow-roll 
parameter 77 as 

n s = l + 2 V , (45) 

where 

, = Sbf^ . (46) 

Imposing the spectral index constraint together with equations pHjl and (jSOJ), we find 
the relations 

m a = -v-, — r> 47 

A m P1 



A 3 = ^ + / A" 3 / 2 — e-^% 48 

m l-exp{-^ c }m 2 

<7 e = 1 . ^4y J 

# rjy/6TrXP n mpi 

With the numerical result for the Euclidean action (|2fij) and the power spectrum nor- 
malization pijl. this gives 

42 . , , rl , m 9.3 x 10 5 , AT ., m 
S E = 7 = (1 - ex V {- V N e }) = — (1 - exp{- V N c }) , (50) 

and, when additionally imposing N e = 60, 77 = —0.025, 

1.3 x 10 8 m 

Se = 7= • (51) 

gV Ampi 

The consistency condition (J24j) with w = w r for our approximation is fulfilled when 

m < 9.7 x 10~ 7 -^=m P1 . (52) 
V A 

Again, we have found that tunneling is preferred for large couplings A and g and small 
values for the mass parameter m, where the small ratio to the Planck scale is imposed by the 
small amplitude of density perturbations. Comparison of the Euclidean actions for the blue 
model ()40j) and the red (JSTj) with iV e = 60 shows that both differ only by a proportionality 
factor which is irrelevant with respect to the level of our approximation. The figure for the 
blue model is therefore almost indistinguishable for the eye when compared to Fig. |21 for 
the red model, which is why it is omitted here. 

We note that for the above models, one should bear in mind that in order to obtain the 
effective potentials (JMj) and ()43j) in the parametric range which allows for tunneling decay, 
one has to require a more than substantial amount of tuning. For hybrid inflation at the 
Electroweak scale, this is discussed in |T4^ I15j. The one-loop correction to the tree-level 
hybrid potential (jSJ is given by 

t/ 1 ( 2 2 2\2, g 2 a 2 -m 2 

^i-ioop = ^ [m -go) log — 2 , (53) 
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where A is an ultraviolet cutoff. Suppose now, we choose the renormalizable counterterms, 
which are the terms up to fourth order in a, in such a way that for o = cr e , we have the 
desired values for m a and A, while the cubic and quartic self-couplings cancel to zero. When 
now expanding the potential around a e down to values of a > cr c , the nonrenormalizable fifth 
order term, which we did not eliminate, induces an additional mass for the field a of order 
gm. This is to be compared with m a as in (jHZj) or (}4?j) . In order for the quantum correction 
to be subdominant, one therefore has to require g <C m 3 /mp l (blue) or g <C m 2 /mp l 
(red), respectively. Comparing with Fig. it is easy to see that tunneling does not play 
any role within hybrid models which do not require the fine-tuning of nonrenormalizable 
operators. This reasoning already strongly indicates that the supersymmetric scenarios 
with radiatively lifted flat directions do not suffer from tunneling either, as we shall see 
explicitly in the next section, number 

The above study shows that inflation exit via tunneling is mostly relevant for inflation 
models with an energy scale below the GUT scale and especially for the intriguing case 
of models where inflation is connected to electroweak physics and therefore within exper- 
imental reach. In fact, a hybrid model has been suggested where the role of the waterfall 
field is played by the Standard Model Higgs field, such that the field content only needs 
to be extended by the inflaton singlet a Since these models apparently also bear 

the potential for successful baryogenesis, the enormous fine tuning of the potential may be 
considered worth the price for a minimal field content. However, our analysis shows that 
for the desired parameters m/\f\ = 246 GeV, and the couplings g, X = O(l) in order to 
allow for strong pre- or reheating, a rapid decay via tunneling is inevitable, such that the 
electroweak hybrid models are not viable, even if fine tuned. 

Another motivation for contemplating low scales of inflation originates from models 
explaining the hierarchy between the Planck scale and a fundamental unified scale by the 
presence of large extra dimensions. Inflation then has to take place at energies below the 
unified scale. Kaloper and Linde [TH] point out that hybrid inflation is a possible way to 
keep the vacuum expectation value of the scalar fields in the effective four-dimensional 
potential for inflation below the fundamental scale. While they agree with Lyth on the 
view that when inflation occurs at the Electroweak scale, one faces a severe fine tuning 
problem, they conclude that models with M ~ 10 11 GeV, g 2 = A = 0.1 fit perfectly in the 
hybrid scenario, somewhat in contradiction with our estimate g -C m 3 /m Pl . Note also our 
comment on the tunneling rate for this parametric range at the end of section 13.11 

3.3 SUSY Hybrid Inflation 

From the general arguments on the irrelevance of the tunneling rate within our approxi- 
mation for models without fine-tuning of nonrenormalizable operators, it is already clear 
that tunneling does not play a role in SUSY- hybrid inflation [TTH |2*U] . These models are 
however of special interest since they rely on rather minimal assumptions and in their sim- 
plest version depend only on a single parameter k [21] , which can be determined [22| from 
the latest observational data .11 3 . They furthermore bear the potential of a successful 

3 Before the WMAP3 data became available, only an upper bound on n could be given. 
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embedding of the Minimal Supersymmetric Standard Model in an inflationary scenario, 
possibly linked to a Grand Unified Theory [201 E31 EH ESJ • Due to the importance of these 
models, we derive here an expression for the Euclidean bounce action Se, although it will 
be large and prohibit tunneling. 

F-term SUSY hybrid inflation is implemented by the superpotential 

W = kS (GG - M 2 ) , (54) 

which leads to the tree-level scalar potential 

V = k 2 \GG -M 2 \ 2 + k 2 \SG\ 2 + k 2 \SG\ 2 . (55) 

The involved fields are complex, where S is a singlet, G a gauged multiplet of dimension 
N and G its conjugate. Vanishing of the D-terms relates the vacuum expectation values 
(G) = (G)*. 

We choose the phase of S to be zero and identify 

a = V2Re[S] . (56) 
Likewise, by a unitary gauge choice, such that (Re[G*]) = we can set 

= v^RefG''] . (57) 
In terms of these fields, the potential (|55jl reads 

V = k 2 M 4 - k 2 M 2 <P 2 + -^k 2 ^ + i«W . (58) 

This is a special case of the more general potential J5J with the replacements 

m 2 = 2k 2 M 2 , 
9 = k, 

X = K 2 . (59) 
The lifting of the flat direction is then induced by the Coleman- Weinberg potential [201 EHj 

* = as?" ((y +M ) ""r^-ip— ) + U _M ) ^r^-ip— 

-r 4lo <43}- < 6 °) 

We consider again the situation where a is close to the critical point, such that we can 
approximate 

5^ « v / 21og2^ K 4 M 3 + 0(M 2 Aa) , (61) 
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where the critical point is at a = a c = V2M and Aa = a — a c . Within this approximation, 
the number of e-folds is 

64vr 3 MAo . . 

N * = -r, ■ 62 

v /21og2«Ww| 1 

Imposing the normalization of the power spectrum ([30)1. we get a relation between k and 

such that we can derive 

Aa = (log 2) 1 / 3 2- 3 / 2 3- 1 / 3 vr- 4 / V/ 3 ( MJ\f 1/3 N e . (64) 
Using (|26|) . we find for the Euclidean tunneling action 

S E = 6.1 (kV^) ~ 2/3 A/" 1/3 iV c = 4800/t~ 2/3 AT 1 / 3 AT e . (65) 

Tunneling therefore does not occur within F-term SUSY-hybrid inflation. 

We have also performed a corresponding study for the D-term model [2U l2Zj ; which 
is more involved due to the additional parametric dependence on the gauge coupling con- 
stant. However, as one can already anticipate from the general arguments about radiative 
corrections and tunneling given at the end of section E21 we find that tunneling is also very 
suppressed in these scenarios. We therefore omit a detailed presentation of the derivation 
of this negative result. 

4 Conclusions 

Imposing the normalization of the scalar perturbation spectrum (f3*D|). it is possible to 
estimate for generic models of hybrid inflation the range of parameters where tunneling 
dominates over the slow-roll evolution of the inflaton fields. In order to calculate the 
Euclidean action Se, we have assumed that the bounce solution follows a straight trajectory 
in the field space spanned by the inflaton and the waterfall field. We have numerically 
obtained the action for a particular set of parameters and then used its scaling properties 
in order to calculate the tunneling rates in parametric regions of small couplings and small 
field values, which are difficult to access numerically. This result is expressed in Eq. 126)1 . 
which we have used to derive constraints on hybrid inflation, arising from the requirement 
that tunneling should be suppressed. The consistency of our approach is verified by a 
comparison with the numerically determined results for the bounce action along the curved 
extremal path in two-dimensional field space. 

Our results are best summarized by the formulas (|40|). (|51|) and by Fig. El Tunneling 
may play a role for models with a mass below 10 12 GeV, but can effectively be suppressed 
by small values of the inflaton- waterfall coupling g and the waterfall self coupling A, which 
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in turn imply large expectation values of the inflaton field during inflation or the waterfall 
field after its end. 

Provided one does not allow for the fine-tuning of nonrenormalizable operators, tun- 
neling never constitutes a problem. In particular, one cannot derive any tunneling bounds 
on the parameters of F- or D-term SUSY models. In contrast, models of electroweak hy- 
brid inflation, which need coupling constants of order one for a sufficient reheating of the 
Universe but require fine-tuning, are completely ruled out since the inflaton would rapidly 
decay through bubble nucleation. Leaving alone the issue of stability of the inflaton po- 
tential with respect to radiative corrections, tunneling decay prohibits the realization of 
hybrid inflation when the vacuum energy and all field expectation values are required to 
be at scales below 10 12 GeV. 
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